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On  the  Error  in  Pade  Approximations  for  Functions 
Defined  by  Stieltjes  Integrals* 


Yudell  L.  Luke 
University  of  Missouri 
Kansas  City,  Missouri 


I.  Summary  and  Introduction 


In  a previous  paper,  the  author  derived  an  infinite  series 
representation  for  the  error  in  the  Gaussian  quadrature  of 


= /> 


(t)f(t,z)dt  where  w(t)  is  the  weight  function  and 


f(t,z)  is  known.  We  suppose  that  f(t,z)  can  be  represented 
by  an  expansion  in  series  of  orthogonal  polynomials  (qn(t)} 
which  is  uniformly  convergent  in  [a,b]  . Thus,  f(t,z) 

• -i  rb 

= 2 c q ft)  , c = h I w(t)q  ft)  f (t , z)dt  . Then  the  error 

n=0  n ln v n n Ja  ui v 

can  be  expressed  in  the  form  2 c2n+2+s82n+2+s  where  the 

s=0 

' s depend  only  on  the  properties  of  qn(t)  • Usually  only 

one  or  two  terms  of  the  series  are  needed  to  achieve  a realistic 

estimate  of  the  error.  Here  asymptotic  estimates  for  c for 

n 

large  n are  employed.  In  the  present  paper,  the  error  analysis 
analysis  is  applied  to  the  situation  where  f(t,z)  = (l+t/z)~*  . 
In  this  instance,  Gaussian  quadrature  of  I and  integrals  simply 


*This  research  was  sponsored  by  the  Air  Force  Office  of 
Scientific  Research  under  grant  73-2520. 
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related  to  I lead  to  approximations  in  z which  occupy  the 

(n,n+r-l)  positions  of  the  Pade  matrix  table  for  I . Three 

important  examples  are  treated.  Two  examples  are  generalized  to 
the  case  where  f(t,z)  = (1+t/z) ~b  . 

II.  Gaussian  Quadrature 


rb 

,a,b)  = I w 
~'a 


(t)f  (t ,z)dt 


where  we  suppose  that  w(t)  and  f(t,z)  are  integrable  over 
the  path  a to  b . We  consider  approximation  of  (1)  by  use 
of  the  Gaussian  quadrature  formula.  To  this  end,  let  {qn(t)} 
be  a system  of  orthogonal  polynomials  of  the  form 


such  that 


I - h i , 


1 if  m = n , and 


6^  = 0 if  m * n . 
mn 


Let  t^  be  defined  by 


= 0 , i 0,1,. ..,n 


and  put 


Ah 
n n 


n+1 ,n+l 


<WWti) 


n a 


3 


Following  Luke  [1],  Gaussian  quadrature  of  (1)  leads  to  the 


formulation 


I(z,a,b)  = In(z,a,b)  + En(z)  , 


In(z,a,b)  = 2 X f (t. ,z) 

J1  i _ n J-  > **  i 


where  E^fz)  is  the  remainder  which  is  zero  if  f(t,z)  is  a 
polynomial  in  t of  degree  < 2n+l  . Under  the  assumption  that 
f(t,z)  can  be  expressed  as  an  expansion  in  the  polynomials 
{ qn ) > which  is  uniformly  convergent  in  [a,b]  , that  is. 


then* 


where 


f(t,z)  = 2 c,qk(t)  , 

k=0  K K 

-i  rb 

= hk  J w(t)qk(t)f (t,z)dt  , 


EnCz^  " sfQ  c2n+2+sg2n+2+s  » 


„(n)  _ hn+l  " An+k+l  _ a2n+2 , 2n+2hn+l 

g2n+2  " aTT  “ A,  2 t12' 

0,0  k=0  k a . 

n+l,n+l 

and  expressions  for  g|n+2+s  » s > 0 * are  detailed  in  the 
source  cited.  There  explicit  forms  for  s = 0,1,2  are  derived 
for  the  classical  orthogonal  polynomials  of  Jacobi,  Laguerre, 
Hermite  and  Bessel, 


L 


*In  references  [1,4],  the  n superscript  notation  on  fi2n+2+s 

was  omitted.  It  should  be  noted,  for  example,  that  the  formula 
for  gk+1  is  not  necessarily  that  for  gk  with  k replaced  by 

k+1  . Here  we  have  added  the  superscript  notation  to  alert  the 
reader  to  this  situation. 
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We  also  have 


n+1 


,2n+2r,.  . 

d fCt,z) 


En(.) 


dt 


2n+2 


t=0 


(2n+2)  a 


a < 0 < b 


n+1 ,n+l 


(13) 


and  clearly  En(z)  5 0 if  f(t,z)  is  a polynomial  in  t o£ 
degree  £ 2n+l  . 

III.  Pade  Approximations 


Let 

F(z)  = 2 v,zk  . 

k=0  K 


(14) 


Let  Pm (z ) and  Qn(z)  be  polynomials  in  z of  degree  m and 
n , respectively,  such  that 

Qn(z)F(z)  - Pm(z)  = 0(zm+n+1)  . (15) 


Then  Pm(z)/Qn(z)  is  the  Pade  approximation  to  F(z)  which 
occupies  the  (n,m)  position  in  the  Pade  matrix  table.  If 
m = n , we  have  the  main  diagonal  entries  and  if  m = n-1  , we 
have  the  first  subdiagonal  entries.  For  further  information  on 
Pade  approximations  with  references  and  numerous  examples,  see 
[2,3]?  See  also  references  [6-9]  ahead. 


*In  12],  the  first  subdiagonal  approximants  are  said  to  occupy 
the  (n-l,n)  positions  in  the  Pade  matrix  table  in  contradis- 
tinction to  the  present  notation  (n,n-l)  . Also  in  [3 ,pp . 493 , 494 ] 
there  are  a few  (obvious)  typographical  errors. 
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IV.  Application  to  Stieltjes  Integrals 


In  this  section,  we  apply  the  results  of  Sections  2 and  3 
to  the  siutation  where 


f(t,z)  * (1  + t/z)'  . 


In  this  event, 

TX  ^ "l 

t r~  „ v*  _ v i,n  _ polynomial  in  z of  degree  n 
n^z,a,D'  h i + 1 . / z rr  ’ 

i = 0 i'  polynomial  in  z of  degree  *n+l 


En(z)  - 


n+1  _ „ , -2n-2-v  „ , 

~2  . 2n+3  a < ? < b • (18) 

an+l,n+llz 


It  follows  that  In(z,a,b)  is  the  first  subdiagonal  Pade 
approximation  to  I(z,a,b)  inhere  f(t,z)  is  given  by  (16). 
Notice  that  if  a and  b are  real  and  z is  positive,  then 
En(z)  is  also  positive. 

We  now  show  how  to  get  the  main  diagonal  Pade  approximation. 


From  (16) , 


and  so 


I*(z 


f(t,z)  = 1 - |(1  + I-)"1  , 


fb  -1 

,a,b)  = I w(t)dt  - z I*(z,a,b)  , 
a 

rb  -i 

(z ,a,b)  = / tw(t)(l  + t/z)  dt  . 


We  apply  our  basic  quadrature  procedure  to  I*(z,a,b)  . In 
effect,  this  is  a generalization  of  Gaussian  quadrature  for  (1) 
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and  (16)  with  an  additional  abscissa  t = 0 , which  goes  by  the 
name  Radau.  For  a discussion  of  generalized  Gaussian  quadrature 
with  arbitrary  additional  abscissae,  see  Luke,  Ting  and  Kemp  [2). 
For  quadrature  of  (20),  let 

P^Ct)  “ " K tk  , (21) 


tw(t)p_  (t)p  ft)dt  = j 6 , 

yrnv  J -’n  mn  ’ 


Pn ( ) = 0 , i = 0,1, ... ,n  , 


F i 

v _ n n_  n _ "n+1, n+1 

ri.n  P^rTTp^TTT  ’ n ' t>nn 


Then 


I(z,a,b)  = I*(z,a,b)  + VR(z)  , 


I*(z,a,b)  = [b  w 


1 n y • 
(t)dt  - i 2 


Z i=0  1+Ui7z 


_ polynomial  in  z of  degree  n+1 


polynomial  in  z~  of  degree  n+1 


where  Vn(z)  is  the  remainder  and 

y (-»•*  - _ ^n+1  . -2n-3v 

nt  3 ^ ;;+A, 2n+3  - °(z  ) » 

bn+l,n+l<z+*) 


a < <(>  < b 


Thus,  I*(z,a,b)  is  the  main  diagonal  Fade  approximation  to 
^j1(z»a»t>)  . If  a and  b are  real  and  z is  positive,  then 


/ 


Vn(z)  < I(z,a,b)  < En(z) 


Since 


(l+t/z)"1  - . z (»r(t/z)k  + C-)rCt/Z.)r,Ci-t/z.)-1  ? . (29) 

k=0 


k -k 


rbtk 

I t w 
->a 


Ct)dt  + {(-) 


r/zr>J^Vw| 


I (z ,a,b)  = Z (-)  Z i t wCt)dt  + {(-)  /z  }/  t w(t)(l+t/z)  dt, 
k=°  *^a  (30) 

Gaussian  quadrature  of  the  integral  in  (30)  leads  to  an  approxi- 


mation which  occupies  the  (n,n+r-l)  position  in  the  Pade  matrix 
table. 

That  Gaussian  quadrature  of  Stieltjes  integrals  leads  to 
Pade  approximations  is  rather  well  known.  Indeed  it  goes  back 
to  Stieltjes  [5] . In  recent  times.  Pade  approximations  of 
Stieltjes  integrals  via  Gaussian  quadrature  together 
with  error  analyses  have  been  considered  by  a number  of  writers. 

We  make  no  attempt  to  give  a complete  bibliography,  but  see  • 

Allen  et  al  [6,7],  Chui  [8],  Karlsson  and  Sydow  [9]  and  the  refer- 
ences noted  in  these  sources.  There  is  an  extensive  literature 
on  Gaussian  quadrature  for  general  integrals  which  we  do  not  quote. 
In  this  connection,  see  Davis  and  Rabinowitz  p.0]  and  Donaldson 
and  Elliott  [jll]  . The  latter  is  especially  valuable  as  it  de- 
velops a unified  approach  to  quadrature  formulas  with  the  error 
expressed  as  a contour  integral  from  which  asymptotic  estimates 
of  the  error  can  be  deduced.  Our  approach  differs  from  all 
other  authors  in  that  we  study  the  error  in  Pade  approximations 
by  use  of  (11).  We  now  turn  to  three  examples. 

V.  Examples 

We  treat  three  cases  for  which  considerable  information  is 
known  about  their  Pade  approximations.  In  what  follows,  we 
make  free  use  of  notations  and  ideas  found  in  my  books  [2,3]|  . 
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In  each  situation  we  treat  the  first  subdiagonal  Pade  approxi- 
mation only. 

Example  1.  Consider 

P n a-c-i/zi  - r(c)  f1  t^Cl-t)^8"1^ 
2rj(l»a,c  1/zj  r (a) r (c-a)  I (1+t/z) dt  » 

J f\ 


R(c)  > R(a)  > 0,  |arg(l+l/z)  | < ir. 

We  have  [2,  v.2,  p.3l], 

oo 

(l+t/z)"1  = I c (z)R^a>3) (t)  , a = c-a-1  , 8 = a-1  , (32) 

n=0  n n 

n«r-ln  /n+B+l,n+l  \ 

c (z)  = n-;-^ ?F  -1/z  , X = ct+8  + 1 = c-1  , (33) 

z (n+X)n  2 1\  2n+X+l  / 

where  R^a,3^(t)  is  the  shifted  Jacobi  polynomial.  From 
[2,  v.l,  p.  237]  , 

, ^n03-X„r  -Un+1)  M . -^X-B-3/2  . 

c (z)  = z(mr)  e _ (l__e  J {1  + 0(n  X) } , (34) 

n (1-e"^)'2"3 


e = 2z  + 1 + 2(z2  + z)!s  I 

where  + sign  is  chosen  so  that  |e  < 1.  This  is  possible  for 

all  z except  -1  _<  z < 0.  From  [l]  , 

g(n)  _ r(n+q+2)r(n+B+2)r(n+X+l)  (n+1)  !r(4n+X+4) 

2n+2  {T (2n+X+2) }2r(2n+X+3) (2n+2) ! 

= 2~X(v/2n)h{l  + CHn'1)}  , 

p Cn) 

g2n+3  _ 2(a-8)  (n+1)  (n+2)  (4n+X  + 4)  _ _f  pUl  ofn-1!} 

"W (2n+X  + l)(2n+X  + 3)((rri+^)  (aB)tl+0(n  )I. 

g2n+2 

We  suppose  that  22n+2+r^g2n+2  = h(a,8){l  + 0(n  *)}  uniformly 
for  all  r 1 where  h(q,B)  depends  only  on  a and  8.  For  the 
special  cases  a = -8  , we  have  verified  this  for  r «=  2.  But 


k 


y 


proof  of  the  general  situation  seems  elusive.  From  (11)  and 
(34) -(36),  we  find  with 

Fn(z)  = r(a)r('c'W  En(z)  » 


(37) 


n 


(Z)  = r(X^l)ize~^2nt3:i(Ue~V^  { + 

~ : - - - 2X-3  - -F^U-  ft  11 


r(a+l)r(6+l)2 


(1-e"^) 


x{l  + (a-B)e"?  + A(a,3,0  + 0(n_1)}  , 

- 2£ 

where  A(a,3,£)  depends  only  on  a,  3 and  £ and  is  0(e  ^) . From 
[2,  v.2,  p.173]  (there  put  a = 1,  cr  = 3 + 1,  p - X and  to  con- 
form to  the  present  notation,  replace  n by  n + 1),  we  have 


(38) 


R. 


n+1 


r i __8+l0_  ^ (2n+3+3)  2a  i 

(z)  = r.(*+l).iTZ ? Cl  + e — ]_  {1  + q^"1)}.  (39) 

r(a+l)r(3+l)2za"i 


-1 


Neglecting  0(n  ) terms,  we  have 


F (z)  {1  (a- 3) e + A(at3,Q} 

TT^TiJ  = (4ze ~^)  ^ (l  + e^)a+li 


(40) 


If  z is  sufficiently  large  so  that  e ^ is  sufficiently  small,  then 


F (z) 
nv  J 


R_ 


U + (a-3) e~^  + A(a ,3,0) 

d + (a-3)e-5  + he~2^{ (a- 3) 2-  X}+  0(e'3^)) 


(41) 


viW 

and  so  within  the  limitations  noted,  the  two  different  error 
analyses  give  the  same  result. 

Example  2.  We  consider  the  complementary  incomplete  gamma  function 


T(v.z)  = 


rV^fV1! 

■(I'-vTJo  TTf 


t vdt 


t7TT  • 


m-V,  JQ 

Considerable  information  on  the  main  diagonal  and  first  sub- 
diagonal Padd  approximations  are  known.  See  [2 , 3]  where  asymp- 
totic results  are  given  to  show  that  these  approximations  con- 
verge, although  the  convergence  is  rather  slow.  Formally,  at 


(42) 


least,  we  have  [2,  v.2,  p.18], 


(1+t/z)"1  = E k(z)I,n  (t), 

n=0  n n 


kn(z)  = {r(n+l-v)}  ^ ^ “ n'z*  VUCn+l-v;l-v;z), 

(t)  = {(l-v)n/nl }1F1(-n;l-v;t)  . 

Here  (45)  is  the  generalized  Laguerre  polynomial.  From 
[2,  v . 1 , p.  264]  , 


L(-v)(t)  , Int) 


hv-h  t/2 


h v 
tt  2n 


-{1  + OCn-1)} 


x cos(2(nt)^  + tt  C*2V-%)  + 0(n~'5)}  , (46) 

with  t fixed  and  n «>.  Also  from  [2,  v.2,  p.200]  or  from 
the  work  of  Wimp  [12] , 

^(z)  = (irz)^(n/z)^v"iez^2exp{-2(nz)'5Hl  + Ofn"2)}, 

z fixed,  | arg  z | £ tt  - e , e > 0 , n 00  . (47) 

Thus  under  the  conditions  stated 

Ln  V‘)(t)kn(z)  = n"sO(exp{-2(nz)^})  , (48) 

and  so  (43)  converges,  though  the  convergence  is  rather  slow. 
However,  the  convergence  is  not  uniform  in  0 < t £ °°  as  can  be 
deduced  from  a known  asymptotic  expansion  for  the  generalized 
Laguerre  polynomial  for  large  n which  holds  uniformly  for 
large  t.  See  Wyman  [.13]  . Hence  for  the  example  (42),  our 
method  of  error  analysis  fails. 

Example  3.  We  now  consider  a form  related  to  the  incomplete 


gamma  function, 
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f Y+ 1“  D _x_i  .1 

1F1(l;X+l;-z)  = (T (X  + l) /2ni}  / epp  A A(l+z/p)  Adp  , 

J y - i°o 

Y > 0,R(X)  > - 1 , | arg(l  + z/y)  | < n . (49) 

This  example  is  important  since  it  illustrates  use  of  Gaussian 
quadrature  and  our  technique  of  error  analysis  for  the  inversion 
of  Laplace  transforms.  From  [2  , v.2,  p.7^,  we  have 

00 

(1+z/p)"1  » S v (z)R  (X,p)  , (50) 

n=0  n n 

vn(z)  = {znun(z)/(n+X)n>  , un(z)  = e”z1F1(n+X;  2n+X+l  ;z)  .1  (51) 

VX’p)  = 2F0C"n'n+X;p’1)  » C52) 

The  latter  is  called  a Bessel  polynomial  and  for  properties  of 
same  needed  in  our  present  study,  see  [l , 2 , 3]  . From  Luke  pL4]  , 
we  have 

u (z)  - „-y  t"*U/C2n+A+l){1  . z2(n+l).(n*X) .oCn'3)}.  (53) 

n 2(2n+X+2) (2n+X+l) 


From. [1] , 


g2n+2 


(-)n+1r(n+X+l) (n+1) ! T (4n+X+4) 
{r(2n+X+2)}2r(2n+X+3) 


f ..  n+l-2n- X + 3/2  -X 
(-)  2 n 


{1  + C^n'1)}  , 


(54) 


«£Jj 

^2n+2^ 


4 (n+1) 2 (4n+X+4) 
(2n+X+l) ( 2n+X+3) 


4n{l  + 0(n  *)  } . 


From  (51)  and  (53) , 


v2n+3  , z(2n+X.2)eZ^-X)/C4^X  + 5)  (4n  + X + 7^ 
V2n+2  (4n+X+4) (4n+X+5) 


{1  + 0(z2/n2)}  . 


(55) 


(56) 


We  suppose  that 
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*’2n+2  + r 2n+2  + r 
e2n+2v2n+2 


= £ (1  + 0(n  1) } , 

r,n  v J * 


(-)T  /z(2n+X  + 2)V 
•,n  ‘ r!  V 4n+X+5  / * 


uniformly  in  r,  r 1.  This  is  true  for  r = 1 in  view  of 
(55,56).  Using  a result  in  [lj , we  can  verify  (57)  when  X = 1 
for  r = 1,2.  However,  proof  of  the  general  statement  appears 
elusive.  If  U^(z)  is  the  error  in  the  first  subdiagonal 
Pade  approximation  for  the  left  hand  side  of  (49) , then  'from 
(11),  (51-54)  and  (57),  we  find 

„ _ (-)n+1Trr(X  + l)z2n+2e"Zr(n+X  + l)(n+l)!{l+0(n~1)>,  (! 


24n+2X+3r (n+1+A/2){r(n+3/2+X/2) }2r(n+2+X/2) 


since 


ez(2n+X+2)/(4n+X+5)  z f = l f (59) 

r=0  r’n 

and  (58)  is  precisely  the  result  given  in  [2,  v.2,  p.191], 

where  to  conform  to  our  present  notation, we  must  replace  n 

and  v by  n + 1 and  X respectively,  and  set  a = 1. 

VI  Extension  to  Generalized  Stieltjes  Integrals 

It  is  of  interest  to  indicate  an  extension  of  the  above 

analyses  to  the  case  where 

f(t,z)  = (1  + t/z)'b  (60) 

where  b is  arbitrary  so  long  as  the  integral  (1)  has  meaning, 
in 

Then/the  quadrature  formula  (8),  the  summation  part  holds  save 
that  f(tj.,z)  is  now  computed  from  (60),  since  the  weights  X^  ^ 
are  independent  of  b.  Now  in  the  expression  for  the  error,  see 
(7)  , (11)  , (12)  and  [l],  the  coefficients  Cj,  depend  on  b,  but 


13 


but  the  coefficients  82n+2  + s ^ePen<^  only  on  the  system  of  orthogonal 


polynomials  and  so  are  independent  of  b.  The  expansion  (60)  in 


series  of  shifted  Jacobi  polynomials  and  in  series  of  Bessel  poly- 


nomials after  the  manner  of  the  corresponding  expansions  for  b = 1 


as  in  Examples  1 and  3,  respectively,  are  readily  deduced  from  re- 


sults given  in  the  cited  references.  With  such  data  in  hand,  (38) 


and  (58)  are  easily  generalized. 


For  a generalization  of  Example  1,  we  have 


2Pl(a,b;c;-l/z)  = 

valid  under  the  same  conditions  as  for  (31) . 


vo  ri+t/ziD 


.Let  Fn(z,b)  be  the  error  in  the  Gaussian  approximation  for  (61) 


such  that  when  b = 1 we  get  the  error  in  the  first  subdiagonal 


Padd  approximation,  see  (38).  Then  under  the  same  assumptions 


which  led  to  (38) , we  find 


Fn(z»b)/Fn(z>1)  = {r(b)}"1{2z(z2  + z)‘?5}b"1{l  + 0(n" 1)  } 


Similarly,  for  a generalization  of  Example  3,  we  have 


1F1(b;X  + l;-z)  = {r  (X+l) /2ni}  / ePp^'^l  + z/p)  "bdp  , 

Jy-ioo 


valid  under  the  same  conditions  as  for  (49).  Noav  let  Un(z,b)  be 


the  error  in  the  Gaussian  approximation  for  (63)  such  that  when 


b = 1,  we  obtain  the  error  in  the  first  subdiagonal  Pade  approxi- 


mation, see  (58).  Then  under  the  same  assumptions  which  led  to 
(58) , we  get 


Un(z,b)/Un(z,l)  = {(2n)b_1/r(b)}{l  + 0(n _1) } . 


In  each  example,  the  dependence  of  the  error  on  b is  quite 


weak . 
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